
    gen Aavg_denom2 = Aavg_denom*1000
    local depvar Aavg_denom2
    label var Aavg_denom2 "Average share denomination"
    gen zz = Ataxpayers_s 
    label var zz "Tax payers per capita" 
    egen x = count(1), by(`depvar' zz) 
    sc `depvar' zz [aw=x] if tag_region==1 || lfit `depvar' zz if tag_region==1, title("")  ///
    ytitle("Average share denomination (NOK)") ///
    ylabel(0(500)2500, angle(horizontal) format(%9.0f))  ///
    xlabel(0.25(0.05)0.4 , axis(1) format(%9.2f))  ///
    legend(order(1 "Region" 2 "Fitted values" )) ///
    saving(BMO_figure3,replace)
graph export .\Output\BMO_figure3.pdf, as(pdf) or(landscape) replace 
